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The importance of measuring the b-quark mass at different scales is emphasized. 
The recent next-to-leading order calculation of three jet heavy quark production 
at LEP and its use to measure mi,(m,z) is discussed. 



1 Introduction 

The question of the origin of the masses of quarks and leptons is one of the 
unresolved puzzles in present high energy physics. To answer this question one 
needs to know precisely their value. However, quarks are not free and their 
mass has to be interpreted more like a coupling than an inertial parameter and 
it can run if measured at different scales. Moreover, in the standard model 
(SM) all fcrmion masses come from Yukawa couplings and those also run with 
the energy. To test fermion mass models one has to run masses extracted 
at quite different scales to the same scale and compare them with the same 
"ruler" . This way, for instance, one can check that in some unified models the 
6-quark mass and the r-lepton mass, although different at threshold energies 
they could be equal at the unification scale. For instance, in fig. |] we plot the 
evolution of the b-quark and the r-lepton masses for both the standard model 
and the mininal supersymmetric standard model (MSSM). We see that the 
MSSM behaves much better than the SM since "unification" of masses could 
happen at much higher energies in the MSSM. 

An important part of the running occurs at energies accessible in present 
experiments. However, the running of fermion masses, although predicted 
by quantum field theory, has not been tested experimentally until now. The 
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Figure 1: Running of the b-quark (dark) and r-lepton (gray) masses for the standard model 
(solid) and the MSSM (dashed), b-quark mass running if given for two values of the strong 
couling constant and the mass. 



reason being that for energies y Q 2 much higher than the fermion mass, m q , 
the mass effects become negligible since usually they are suppressed by m 2 /Q 2 . 
But this argument is not always correct. In fact, in the SM, the partial decay 
width of the Higgs-boson into b-quarks is proportional to m 2 (m kJ because the 
Yukawa coupling of the Higgs-boson. In addition, it is not true 111 for jet cross 
sections which depend on a new variable, y c (the jet-resolution parameter that 
defines the jet multiplicity) and which introduces a new scale in the analysis, 
E c = \J Q 2 Vc- Then, for small values of y c there could be contributions coming 
like m 2 /E 2 — (m 2 /Q 2 )/y c which could enhance the mass effect considerably. 
In addition mass effects could also be enhanced by logarithms of the mass. For 
instance, the ratio of the phase space for two massive quarks and a gluon to the 
phase space for three massless particles is 1 + 8(m 2 /Q 2 ) \og(m q /Q). At Q 2 = 
m 2 z and for the bottom quark this gives a 7% effect, for mj, = 5 GeV and a 3% 
effect for m& — 3 GeV. n Thc high precision achieved at LEP makes these effects 
relevant. We showed □ that b-quark mass effects in three-jet production at 
LEP are large enough to be measured and, in factjjjhfiy were already observed 
in the tests of the flavour independence of a s (mz )□'□• Therefore, these effects 
could be used to measure the b-quark mass at LEP. 

The observable proposed—as a means to extract the bottom-quark mass 
from LEP data was the ratio on 

3 " it („ e )/r< ■ (1) 



In this equation T^j(y c )/T q is the three-jet fraction of Z-decays and q denotes 
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the quark flavor. Obviously, the ratio R^ d depends on the jet-clustering algo- 
rithmLl used to define the jets. In this ratio and at the leading order (LO) the 
quark mass effects can be as large as 1% to 6%, depending on the values of 
the mass and the jet-resolution parameter, y c . 

Since the measurement of R^ d is done far away from the threshold of b- 
quark production, it can, in principle, be used to test the running of a quark 
mass as predicted by QCD. However, the leading order calculation does not dis- 
tinguish among the different definitions of the quark masscl: perturbative pole 
mass, Mb, running mass at -Mf,-scale, or running mass at m^-scale. There- 
fore, to distinguish them it is necessary to use a complete next-to-leading 
order (NLO) calculation of threedet^ralios including quark masses which w. 
not available until very recently 0{]€rulij. Here we overview the calculation 
used Jijj-thij DELPHI Collaboration to extract the 6-quark mass at the mz 
scaleBBB 
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2 Jet ratios with heavy quarks at NLO 

The decay width of the Z-boson into three jets with a heavy quark can be 
written as follows 

1% = ^f^f [g v H v (y c ,r b )+g A H A (y c ,r b )] , (2) 

where g is the SU(2) gauge coupling constant, cw and syy are the cosine and 
the sine of the weak mixing angle, gv = — 1 + 4/3s^ and g A = 1 are the vector 
and axial-vector coupling of the Z-boson to the bottom quark and a s is the 
strong coupling constant. Functions Hy^(y c ,rb) contain all the dependences 
on y c and the quark mass, Tb — (Mb/mz) 2 , for the different algorithms. These 
functions are computed perturbatively as an expansion in a s and can also be 
expanded in rb for r& << 1. At the NLO we have contributions to the three-jet 
cross section from three- and four-parton final states. One-loop three-parton 
amplitudes are both infrared (IR) and ultraviolet (UV) divergent. However, the 
UV divergences are removed after renormalization. The four-parton transition 
amplitudes are also IR divergent but the sum of both, four-parton and three- 
parton, j-cpntributions is IR finite. We use the so-called phase space slicing 
methodt-3 to obtain the remaining finite results, that is, the functions Hy(A) 
in eq. (|^) at order a s . These results have beeu-checked in different ways, in 
particular we checked that the massless result clO is nicely recovered when we 
take the limit Mb ^Xl—These results have also been checked independently by 
two different groups nr _. 
Combining eq. (|l|), eq. (||) and using the known expression for F b Eni3 we 
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write R^ d as the following expansion in a s 



R b 3 d = 1 + r b (b + ^b^j , 



(3) 



where the functions bo and b\ are given by an average of the vector and axial- 
vector parts of the Z- widths, weighted by cy = g v /(g v + g\) an d ca = 
g A /(g v + 9 \) respectively. It is important to note that, because of the particu- 
lar normalization we have used in the definition of R^ d , most of the electroweak 
corrections cancel. 

Although intermediate calculations have been performed using the pole 
mass, we can also re-express our results in terms of the running quark mass 
by using the known perturbative expression M 2 = m 2 (/x)[l + 2a s (fi)/ir (4/3 — 
log(m 2 //i 2 ))] . We obtain 



R b 3 d = l + f b (»)[b + 



26olog^f 
A 1 



(4) 



where rb{^) = ml(fi)/rri^ and 

h = h 



b [8/3-21og(r 6 )] 



(5) 



rb(/z) can be expressed in terms of the running mass of the &-quark at fi = mz 
by using the renormalization group. At the order we are working ffe(/i) = 
fb{m z ) (a s (m z )/a s (/i)) _47o//3 ° with ot s (ji) = a s (m z )/(l + a s (m z )/3ot) and 
t = log( A1 2 /m|)/(4 7 r), 0o = ll- 2JVW3, N f = 5 and 7 o - 2. 

At the perturbative level eq. (g) and eq. (|]) are equivalent. However, 
they neglect different higher order terms and lead to different answers. The 
difference gives an estimate of the size of higher order corrections. 

Here we present only the results for the DURHAM algorithm E3, which 
gives smaller radiative corrections and was the one used by the DELPHI Col- 
laboration in its analysis. 

The function bo describes the corrections due to the quark mass at LO 
and has an almost negligible residual mass dependence. A fit, b a = X)r- n ^n™^ 



log™y c , to the complete results for the DURHAM algorithm gives. n, 



=o "-o 

(0) 



-10.521 , k ( 1] = - 
The function 



,4.4352 , k 



(2) 
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-1.6629 . 



gives the NLO massive corrections to R\ . It is im- 
portant to note that b\ contains significant logarithmic corrections depend- 
ing on the quark mass. We take them into account by using the form b\ = 
k[ 0) 3 log(j/ c ) + k$ log(r b ) in the fit. For the DURHAM scheme we obtain: 



k { ° } = 297.92 , k\ L > = 59.358 , k 



(i) 



46.238 . 
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Figure 2: NLO results for R^ d (DURHAM) written in terms of the pole mass, Mi (dashed), 
or the running mass mi,(mz) (solid) at LO (gray) and at the NLO (dark). The dotted line 
gives the NLO result in terms of the running mass at an intermediate scale m;,(10 GeV). 
We use as a starting point m(,(m(,) = 4.13 to obtain both Mf, and m,j,(mz). 

In fig. H we present R^ d as a function of y c for DURHAM. To compute it, 
we use the low energy measurement of the b-quark massllj, m^ms) = 4.13. 
Note that choosing a low value for /x in the NLO predictions written in terms 
of the running mass makes it closer to the LO result written in terms of the 
pole mass, while choosing a large /i makes the result approach to the LO result 
written in terms of the running mass at the mz scale. 

3 ra\y(mz) from LEP data 

Since R^ d has been measured to good accuracy by DELPHI El one can use 
eq. (^) and the relationship between mb(fi) and mb(mz) to extract m(,(mz). 
On the other hand, one could use eq. (||) to obtain directly the pole mass, Mb, 
and then use the relationship between the pole mass and the running mass and 
the renormalization group to obtain also mb{mz)- Both results are slightly 
different and have a residual dependence on the scale /i. The difference gives 
an estimate of the errors due to the unknowledge of higher order corrections. 
In fig. ^ we plot, as a function of fi, the values of mb{mz) obtained from 
R-3exp = 0.971, by using the two methods. The most conservative estimate of 
the theoretical error is obtained by taking the spread of the results at /x = mz- 
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Figure 3: Extracted value of ra^imz) if = 0.971 as a function of the scale jj, by using 

both the pole mass result and the running mass result. 

The final result, after including also statistical errors and the errors due to 
the uncertainties in theJaadronization corrections has been presented by DEL- 
PHI in this conference 113. The obtained mass of the. bottom-quark, rrib(mz), 
measured from the three-jet decay of the Z-bosoni^l is fully compatible with 
the value obtained from low energy determinations 113 after using the renormal- 
ization group. This provides, for the first time, a nice check of the quark mass 
sector of QCD in a very wide range of scales. These results can probably be 
improved by understanding better the scale dependence of the results, by re- 
summation of large logs, by using other observables with a softer dependence 
on the scale, by reducing the hadronization uncertainties and finally by in- 
cluding all LEP data. We think that these studies can provide a rather precise 
value of rrib(mz)- 
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